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Abstract: We are now operating the 500 m 2 Yangbajing air-shower core (YAC-II) array near the center of the Tibet 
air-shower array ( Tibet-Ill ) to observe cosmic-ray chemical composition at the knee energy region since February 201 1 . 
The first step of YAC, called YAC-I, containing 16 detector units, was operated from May, 2009 to February, 2010. In 
this paper, we used the YAC-I and Tibet-Ill coincident data set obtained from May, 2009 through January, 2010 to present 
the electromagnetic spectrum of air shower cores at around 10 15 eV energy region. The effective live time is calculated 
as 100.5 days. We would like to report the comparison of our experimental data with MC model prediction in this paper. 
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1 Introduction 



Direct measurements of the primary cosmic rays (CR) 
with energies higher than 10 15 eV are difficult due to their 
low flux and the limited detector acceptance of the on board 
satellite or balloon experiment. Instead, their properties are 
reconstructed from the measurements of the extensive air 
showers (EAS) they produce in the atmosphere. The recon- 
struction of EAS events is based on Monte Carlo hadronic 
interaction models of the air shower development, which 
are based on the knowledge obtained from the accelera- 
tor hadron-nucleus collision experiments. Since accelera- 
tor experiment can not provide all information that cosmic 
ray studies need, some extrapolation to higher energies and 
to un-reached phase space is inevitable that induces uncer- 
tainty in the explanation of AS phenomenon. 

It is well known that the produced particles in the most 
forward region of hadronic interactions are most responsi- 
ble for the AS development, and the most forward region 
is the dead-corner of conventional collider experiments. It 
is also known that the high energy particles in the AS core 
region are most sensitive to the forward region particle pro- 
ductions. Because of the advantage observing EAS cores in 
the high altitude, a new hybrid experiment was constructed 
and operated in Yangbajing, Tibet. 

In present work, we report the checking of hadronic 
interaction models by observing EAS cores at the energy 
region of 10 15 eV using Yangbajing Air shower Core de- 
tectors (YAC-I) and the air-shower array (Tibet-Ill). 

2 The Tibet hybrid experiment 

The Tibet hybrid experiment consists of air-shower array 
(Tibet-Ill) and YAC-I. The Tibet-Ill array [1] consists of 
733 scintillation detectors (0.5 m 2 each). Fast-timing de- 
tectors are placed with 7.5 m spacing and density detectors 
are placed with 15 m spacing. An event trigger signal is 
issued when any four-fold coincidence occurs in FT coun- 
ters with each of them recording more than 0.6 particles. 
The primary energy of each AS event is determined by the 
air shower size (N e ) which is calculated by fitting the lat- 
eral particle density distribution to the modified Nishimura- 
Kamata-Greisen (NKG) structure function. 

YAC-I that consists of 16 EAS core detectors is shown in 
Fig.l and for the brief description see [2], which has started 
data taking since May, 2009. YAC-I is located near the 
center of the Tibet-Ill air-shower array, operating simulta- 
neously with Tibet-Ill. For the coincident events Tibet-Ill 
provides the total energy and the direction of air showers 
and YAC-I observes high energy electromagnetic particles 
in the core region. 

If any one of YAC-I detectors makes a trigger signal that 
corresponds to at least 20 MIPs' incidence, all ADC data 
from all YAC-I units are recorded. Also the trigger signal is 
sent to DAQ system for AS array. ADC modules of YAC- 
I are calibrated every 4 hours. ADC pedestal values are 
measured every 10 minutes. Each DAQ system has GPS 



Figure 1 : Schematic view of the Tibet-Ill air-shower array 
and YAC-I array. 

clock module independently. The matching between YAC 
data and AS data is made using coincidence of GPS clocks 
and trigger tag to AS array. The coincidence condition of 
GPS is about 1 [is [3]. 

In the following analysis, we present our results based 
on the YAC-I data and AS data. 

3 Simulation and Analysis 

A Monte Carlo simulation has been carried out on 
the development of EAS in the atmosphere and the re- 
sponse in YAC-I. The simulation code CORSIKA (version 
6.204)[4] including QGSJET2 and SIBYLL2.1 hadronic 
interaction models are used to generate AS events. The 
assumed primary cosmic -ray composition in MC is based 
on Non-Linear Acceleration (NLA) model (about details, 
please see [5] [2]). The factional contents of the assumed 
primary cosmic -ray flux are listed in Table 1. Primaries 
isotropically incident at the top of the atmosphere within 
the zenith angles from to 60 degrees are injected into the 
atmosphere. The minimum primary energy of this simula- 
tion is set at 1 TeV. Secondary particles are traced to the al- 
titude of 4300 m till 300 MeV. For each simulated AS event 
that reaches the observational level, its core is dropped ran- 
domly onto an area of 52.84 m x 52.14 m, which includes 
the marginal space of 25 m outside the each side of detec- 
tors. MC simulation shows that the core resolution is better 
than 2 m if taking the N\, weighted center as the AS core. 
The electromagnetic showers in the lead layer induced by 
electrons or photons that hit any detector unit of the array 
are treated by a subroutine that is based on the detector 
simulation code EPICS (version 8.64)[6]. 

Normally, the following quantities of YAC are used 
to characterize an EAS core event: The number of shower 
particles hitting a detector unit is called 'burst size' (Nb). 
When the burse size of a detector unit is higher than 200, 
this unit is defined as a 'fired' one. We also call the total 
burst size of all fired detector units as Nj,, the maximum 
burst size among fired detectors as A^ op . 
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Table 1: Fractions of components in the assumed primary 
cosmic-ray spectrum of the NLA model. 
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Table 2: 


The fraction of the components after the event selection. 
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Figure 2: The distribution of primary energy of the sample 
of two data sets. 



Table 3: The number of events of 2 selected samples. 
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3858 


2888 


317 



We can obtain different event samples that have different 
average primary energy and different sample size by using 
different threshold of N e . Therefore, we can see how some 
physics quantities change with energy simultaneously. We 
obtain two data sets by imposing the following conditions: 

(1) N b > 200, N hit > 6, 7V* op >1500, 10 5 < N e < 
5 x 10 5 ; 

(2) N b > 200, N hu > 6, iV* op >1500, JV e >5x 10 5 ; 
Fig. 2 shows the primary-energy distribution of these two 
data sets. The mode energy as known from the Monte Carlo 
is 260 TeV and 1800 TeV, respectively. 

We sampled 1.8 x 10 10 and 1.17 x 10 10 primaries for the 
QGSJET2 and SIBYLL2.1 model, respectively. The num- 
ber of events of two data sets selected under the above con- 
ditions can be seen from Table 3. The simulated data were 
analyzed in the same manner as in the procedure for the ex- 
perimental data analysis. In present paper, we used the ex- 
perimental data set obtained from May, 2009 through Jan- 
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Figure 3: The comparison of air-shower size N e between 
MC and experimental data. 
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Figure 4: The spectrum of the total burst size ^ 
obatined by MC and experimental data at 260 TeV (a) and 
1800 TeV (b) energy region, respectively. 



uary, 2010. An event coincidence between AS events and 
YAC-I events is made by their arrival time. Deadtime cor- 
rection of 12% for AS trigger system and 15% for YAC- 
I trigger system are taking into account. The data sam- 
ple coming from successful coincidence corresponds a live 
time of 100.5 days. The total number of events which is 
selected under the above two conditions from the experi- 
mental data is also listed in Table 3. 



4 Results and Discussion 

Since our Monte Carlo simulation is started from 
1 TeV, in order to normalize MC data and experimental 
data, we need to know the integral intensity of all particles 
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Figure 5: The flux ratio of the absolute intensities of the to- 
tal burst size Y^, Nb obatined by MC and experimental data 
at 260 TeV (a) and 1800 TeV (b) energy region, respec- 
tively. 



5 Summary 

The shape of the distributions of ^ Nj, is consistent 
between the YAC-I data and simulation data in these two 
cases, indicating that from 260 TeV to 1800 TeV, the 
particle production spectrum of QGSJET2 and SIBYLL2.1 
may correctly reflect the reality within our experimental 
systematic uncertainty of a level about 10%. 

But note that, NLA composition model used a steeper 
He spectrum [5], our results are still affected by the 
composition model used, comparing with the new results 
from PAMELA and CREAM. Enhancing He spectrum 
may change the results. The bending energy of p and He 
spectra may also be an important factor. It is also noticed 
that, seen from Table 2, the 2 data samples have different 
composition. Therefore, at present stage, it is not simple to 
make a conclusion. A further study is needed and is going 
on. 

The above results show that taking the priority of high 
altitude (like Yangbajing) an EAS core event sample can be 
obtained with high statistics by using YAC type detector, 
and the hadronic interaction models can be checked. YAC- 
II has been constructed and start data taking since August 
1st, 201 1, the more results will be expected. 



of cosmic rays at Eo> 1 TeV: starting from Horandal's 
spectra of each composition[7], we improve the major 
8 ones (P, He, C, O, Ne, Mg, Si, and Fe) by the newest 
measurements [8][9][10]. The resultant integral intensity: 
I(>1 TeV) = 0.139 cm- 2 s- 1 sr- 1 with the error +0.0013, 
-0.0012 coming from the error of the index of each the 8 
spectra. 

Fig. 3 is the comparison of air-shower size (N e ) 
between MC and experimental data which are normalized 
by number of events. It shows both MC models produce 
air-shower size distribution consistent with experimental 
data. 

Fig. 4 is the comparison of the total burst size Y Nb 
which are normalized by number of events between MC 
and experimental data at 260 TeV and 1800 TeV energy 
region, respectively. Yl Nb should depend sensitively on 
the inelastic interaction cross section, the inelasticity, and 
particles produced in the forward region. It shows that 
these two hadronic interaction models have the same shape 
with experimental results. 

Fig. 5 is the flux ratio of the absolute intensities between 
MC and experimental data. It shows that both QGSJET2 
and SIBYLL2.1 give about 40% lower flux. 
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